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Abstract 



We study the effects of dynamical quarks on the static quark potential at 
distances shorter than those where string breaking is expected. Quenched 
calculations and calculations with three flavors of dynamical quarks are done 
on sets of lattices with the lattice spacings matched within about one percent. 
The effect of the sea quarks on the shape of the potential is clearly visible. We 
investigate the consequences of these effects in a very crude model, namely 
solving Schroedinger's equation in the resulting potential. 
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I. INTRODUCTION 



One of the conceptually simplest things that can be computed in lattice QCD is the 
potential between a static, infinitely heavy quark and antiquark [I]. This potential is used 
in phenomenological approaches to hadron physics, and is also commonly used in lattice 
studies to determine the lattice spacing, or overall energy scale. In this work, we explore 
the differences between the potential computed in the quenched approximation and in full 
QCD. 

In general, there are two main effects. The more dramatic effect of dynamical quarks 
is the breaking of the string at large distances caused by shielding of the static sources 
by dynamical quarks. This string breaking is difficult to see in straightforward potential 
calculations, although studies using both "string" type operators and "meson-meson" type 
operators show the expected behavior 0. We do not consider it here. The second effect 
is that at shorter distances the dynamical quarks will change the shape of the potential. 
Roughly speaking, this happens because the quarks reduce the QCD beta function. This 
means that as we go to shorter distances, the coupling constant in full QCD will not decrease 
as quickly as in quenched QCD, so the force between the static sources will be larger in full 
QCD. This effect is the subject of this paper. 

This difference is responsible for a shift of the hyperfine splitting in quarkonia. In poten- 
tial models, the splitting is proportional to the strong coupling constant and to the square 
of the wave function at the origin, which is greater in the full theory because the coefficient 
of the beta function is smaller ||. The same effect has been discussed as an explanation 
for differences in heavy-light pseudoscalar decay constants, fs and /b s , between full and 
quenched QCD [|J. The idea is that the deeper full QCD potential at the origin will increase 
the wave function of the light quark in a potential model, hence increasing the amplitude 
for it to annihilate with the heavy quark. 

When one describes the effects of the number of flavors on some physical quantity, one 
is making an implicit assumption that some other physical parameter is being held constant 
as the number of dynamical flavors is varied. We are also interested in a variation of this 
problem — we would like to use a quantity determined from the heavy quark potential as 
the parameter which is held constant while other physical parameters are compared. To do 
this, we have to keep in mind that quenched QCD and full QCD with different numbers of 
dynamical quark flavors are different theories. The potentials of these theories have different 
shapes. It is impossible to determine the relative length scales in the different theories by 
matching the entire potential. Instead, one must make a choice — for example, determining 
the lattice spacing from the string tension a or from Sommer's parameter ro @. If the 
couplings of full and quenched simulations are tuned by matching one such parameter, they 
will differ in the other parameters. 

In order to isolate the effects of dynamical quarks, we have done a series of simulations 
of quenched QCD and three-flavor QCD with a range of quark masses, with the gauge 
couplings tuned to match the lattice spacing, using one of the possible choices of definitions 
of the length scale. We used a Symanzik improved gauge action and, in the full QCD runs, 
an improved Kogut-Susskind quark action that greatly improves the scaling of physical 
quantities p|-|l3"||. The full QCD simulations were done with an optimized version of the 
updating algorithm for general quark actions described in Ref. 0. Most of the results 
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presented here come from a series of runs with three degenerate quarks with masses ranging 
from approximately the strange quark mass to eight times this mass. We also include some 
results from shorter runs with two lighter flavors and one quark mass fixed at the strange 
quark mass. Some of the parameters of these simulations are tabulated in Table |. Lattice 
spacings are determined with fractional errors ranging from 0.2 percent for the quenched 
run to 0.6 percent, and the lattice spacings are matched with a fractional RMS spread of 
0.6 percent. 
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TABLES 



TABLE I. Simulations used in this paper. The first six columns are the sea quark mass(es), 
the gauge coupling, the number of simulation time units in the run, the number of configurations 
for which the potential was calculated, the size of the molecular dynamics time step, and the 
approximate number of conjugate gradient iterations used in the computation of the fermion force. 
The remaining columns give r\ , tq and a in units of the lattice spacing, and the lattice spacing in 
fermi, using T\ = 0.36 fm to set the scale. The quenched runs were done with a combination of 
overrelaxation and heat bath steps, using ten sweeps between samples, with each sweep consisting 
of four overrelaxation steps and one heatbath step. 
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The static quark potential including effects of sea quarks has been calculated by the 

1|, UKQCD and CPPACS collaborations and by 



SESAM/TxL 0, SCRI 
Tamhankar and Gottlieb 
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The SESAM calculation used the plaquette gauge action 
and two flavors of Wilson quarks, with a fixed gauge coupling (6/g 2 = 5.6) for the dynami- 
cal runs. The SCRI group used two flavors, with both Wilson and Kogut-Susskind quarks, 
at fixed gauge couplings. The UKQCD calculation used the plaquette gauge action and two 
flavors of quarks with the Wilson-clover action and the gauge coupling was tuned to match 
the lattice spacings. The CPPACS calculation used an improved gauge action with two 
flavors of Wilson-clover quarks, varying the sea quark mass at fixed gauge coupling. Both 
the UKQCD and CP-PACS groups observed an increase in the coefficient of the Coulomb 
term due to the dynamical quarks. Because we are using three flavors of sea quarks, our 
results cannot be quantitatively compared with these two flavor results. 



II. THE SIMULATIONS 



As mentioned above, a choice must be made when deciding how to tune the lattice 
spacing. The most common choice in recent work is to use the parameter tq defined by 
r 2 F(r ) = 1.65, with r ~ 0.50 fm. This choice was motivated by the observation that 
phenomenological studies of the potentials give consistent results for the force at this distance 
||. However, as noted in the original work [[|, in principle r 2 F{r) can be set to any number. 
In these simulations, we are faced with the challenge of tuning the coupling with minimal 
effort, so it is important that the length scale we choose be accurately measurable in a small 
simulation. This motivates us to choose a variant which we unimaginatively call r 1; defined 
by r\F(ri) = 1. Figure [I] shows the fractional accuracy with which the lattice spacing could 
be determined in a short simulation for different choices of the constant in r 2 F(r) = C. 
These are plotted as a function of the distance at which the condition is satisfied. The 
square in this plot indicates the conventional choice C = 1.65, and the octagon shows our 
choice C = 1.0. The arrow on the right side shows the accuracy obtained from the string 
tension. The reason for the improvement in accuracy is that the potential is more accurately 
determined at shorter distances. However, at very short distances the potential is dominated 
by the —a(r)/r part, and in this regime a length scale could only be determined by the small 
dependence of a(r) on r. Thus, the optimal choice for C satisfies r 2 F(r) = C at roughly 
the distance where the force turns over from Coulomb-like to or. 
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FIGURES 
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FIG. 1. Accuracy of a lattice spacing determination from the static potential as a function of 
the distance at which the force is determined. The square is the point where r 2 F(r) = 1.65 and the 
octagon the point where r 2 F(r) = 1. The arrow is the accuracy of a string tension determination. 

We computed the static potential by gauge fixing the lattices to the Coulomb gauge and 
then just computing products of the temporal links. This is effectively a smeared spatial 
string operator, smeared over the distance scale at which the QCD dynamics correlates the 
links, but we have made no further attempt to optimize our static pair creation operator. 
The potential was determined from the ratios of "loops" at time separations four and five 
(five and six for the a ~ 0.10 quenched run). Ratios of distances three and four gave results 
which in some cases differed from the ratios at times 4-5 by more than the statistical error, 
while inclusion of larger distances in a fit did not appreciably improve the statistical errors. 
After the potential was computed at all spatial separations, separations related by lattice 
symmetries were combined. Then the potential was fit to the form V(r) = C — a/r + or. 
This form gave good results for separations from (1,1,0) to the largest values for which 
we could compute the potential. Our quoted results are obtained to fits from all points at 
spatial distances from \/2a to 6a, with the exception of the a « 0.10 fm run, where distances 
up to 7a were used. r and r\ were derived from a and a. Such a simple approach to fitting 
the potential is only possible because the improved action suppresses most of the lattice 
artifacts. Figure [| shows the potential for the quenched run, where our statistics are largest, 
together with the fit that was found. A remaining lattice artifact is visible at distance two 
(separation (2, 0, 0)), where the potential is visibly below the smooth curve. An elegant test 
for lattice artifacts is to compare the two points at distance three, at separations (3, 0, 0) 
and (2,2,1), which are shown in an inset in Fig. [|. For the quenched potential there is 
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a small discrepancy here, with the "on-axis" point lower. For most of our full QCD runs 
this discrepancy is smaller than the statistical errors. For the quenched case we have done 
a shorter simulation at a smaller lattice spacing, a pa 0.10 fm to check for dependence on 
the lattice spacing. This run was done on 28 3 x 96 lattices at 10/g 2 = 8.4. Figure § shows 
the short distance part of the potential for the quenched theory at the two lattice spacings, 
showing excellent agreement. Since the main point of this work is to isolate the effects 
of dynamical quarks, we have not included a systematic error from the remaining lattice 
artifacts in our quoted error estimates. First, we think that this error is small. But more 
importantly, we are comparing quenched and full QCD lattices at matched lattice spacings, 
and these systematic errors will largely cancel in the differences of the two. (All of our runs 
show essentially the same discrepancy at distance two, for example.) 

As detailed in Table |, we were able to choose values of the gauge coupling to match t\ in 
the different runs with a fractional RMS spread of 0.6 percent, with the final lattice spacings 
known with accuracies ranging from 0.2 percent for the quenched run to 0.6 percent for the 
lowest mass run. 

Our procedure for determining the lattice spacing differs from that used in Ref. |J (see 
also [|1^]), where an interpolation is done using values of the potential along one axis, using 
a tree level calculation of the potential to correct for lattice artifacts. We have compared the 
results of the different methods using our large 10/g 2 = 8.0 quenched sample. In Table |Tl 
we tabulate r /a and r 1 /a calculated using our fitting procedure and by the procedure of 
Ref. [|J applied along the lattice directions (1,0,0), (1,1,0) and (1,1,1). The results are 
consistent, with r\/a better determined than r^/a for both methods. The fitting method 
produces smaller statistical errors simply because it makes use of the potential at all possible 
separations. 

All of the errors quoted in this paper are statistical errors obtained with the jackknife 
method. 



TABLE II. Accuracy and consistency of the lattice spacing determination from fitting all 
points and from the method of Ref. || applied along three different axes. For the method of 
Ref. H the second column shows the direction along which interpolation was done. 
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FIG. 2. The static quark potential in quenched QCD, with 10/ gf mp = 8.0. Octagons are from 
time separation four to five, and diamonds from time separation three to four. The solid line is the 
Coulomb plus linear fit to the potential. Insets show the remaining lattice artifacts, as explained 
in the text. A constant has been subtracted from the potential to set V(r\) = 0. The vertical lines 
are at r\ and ro- 
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1 2 

t/t 1 

FIG. 3. The short distance part of the static quark potential in quenched QCD for two different 
lattice spacings. Octagons are from the W/g 2 = 8.0 run, and diamonds from the 10/ g = 8.4 run. 
The solid lines are the Coulomb plus linear fit to the potentials, plotted in units of r\. The upper 
two sets of vertical lines show the lattice spacings — a, 2a ... — for the two runs, and the lower 
set is at spacings of 0.1 fm, where the overall length scale was set by tq = 0.50 fm. 
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III. EFFECTS OF SEA QUARKS 



In Fig. |], we plot the quenched potential and the potential with three degenerate flavors 
of quarks at am q = 0.05. This is approximately the mass of the strange quark. Both the 
distance scale and the potential are plotted in units of r 1; and a constant has been subtracted 
from the potential so that it is zero at r±. Since r\ was determined from this potential, this 
results in matching the slope of the potentials at r = r 1; so the fits are tangent at this 
point. Away from the point where the force is matched, the potentials have different slopes, 
as expected. This graph contains three "rulers" . The upper sets of vertical lines show the 
lattice spacings for the dynamical and quenched runs, and the lower set is at spacings of 0.1 
fm, where the overall length scale was set by ro = 0.50 fm in the quenched run. 
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FIG. 4. The static quark potential for quenched (octagons) and three flavor (diamonds) QCD, 
in units of r\. The solid lines are fits to the Coulomb plus constant plus linear form. The lattice 
spacing was matched using r\ as described in the text. As in Fig. ||, the upper two rulers show the 
lattice spacing in the two runs, and the lower one shows units of 0.1 fm. The inset expands the 
area shown by the box. 
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As a quantitative measure of the change in the shape of the potential due to the dynamical 
quarks, we may plot dimensionless quantities such as r y/a or r 1A /cx. Figures [| and || show 
these two quantities as a function of quark mass, where {m^/rrip) 2 is used as an indicator of 
the quark mass. The effect of the sea quarks is clear, and it should be noted that the effect 
is appreciable even for fairly heavy sea quarks. (The dependence of these quantities on the 
sea quark mass emphasizes that only one dimensionful quantity can be matched in tuning 
full and quenched simulations to match the lattice spacing.) Figure [5] can be compared with 
Fig. 1 of Ref. flUf ) which plots the same quantity with two flavors of Wilson quarks. We note 
in passing that an extrapolation of Fig. [5] to light quark mass is consistent with commonly 
used values for r and y/a: 0.50 fm x 440 MeV = 1.12. 
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FIG. 5. Effects of dynamical quarks on the shape of the potential. Here we plot rQy/a as a 
function of the quark mass. The quenched points are at the right, with the octagon coming from 
the W/g 2 = 8.0 run and the cross from the 10/<? 2 = 8.4 run. The remaining octagons are full QCD 
runs with three degenerate flavors, and the squares are full QCD runs with two light flavors and 
one heavy. 
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FIG. 6. The same as Fig. ^j, except we plot v\yfa. Physically, the difference is that this is 
sensitive to shorter distances than the previous plot. 



IV. PHENOMENOLOGICAL CONSEQUENCES 

The larger quark- ant iquark force at small distances in full QCD has been suggested as 
a qualitative explanation for some of the differences between quenched and full QCD that 
have been found in simulations. For example, there is some indication that dynamical quarks 
increase the value of fs and fs s fl20|- j23|l , which can be understood as a larger wave function 
at the origin due to the deeper short distance potential well in full QCD [|]. Also, dynamical 
quarks are expected to increase the spin splittings in heavy quarkonia || 24,|!J, which can 
again be considered to result from a smaller wave function. 

To make a crude estimate of these effects, we have solved the nonrelativistic Schroedinger 
equation for the ground state in the static potentials. Of course, we don't advocate this 
as a replacement for a real solution to QCD, but in the spirit of trying to understand the 
differences between quenched and full QCD, it is a worthwhile exercise. Specifically, we used 
the fitted Coulomb plus linear plus constant form for the potential. Alternative approaches, 
where we connected the points in the potential by straight lines or solved a discrete lattice 
version of the Schroedinger equation, gave similar results. We show the "Coulomb plus 
linear" results because we expect that this is closer to continuum physics than the discrete 
lattice potential, which is finite even at r = 0. However, it must be remembered that this 
1/r term is phenomenological, not fundamental. (We are working at distances too large 
to identify it with a perturbative potential, and too short to identify it with a universal 
correction to the linear potential.) For light valence quark masses the wave function extends 
beyond the region where we have fit the potential and into the region where we would expect 
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string breaking to show up. By using the fitted form of the potential, we are ignoring this 
effect. Conversely, for heavy valence quark masses the wave function is small, and sensitive 
to the treatment of the short distance potential. However, with our matched quenched and 
full QCD lattices, we expect this effect to be the same for the quenched and full potentials. 
Note that we can use our small sample of quenched lattices with a ~ 0.1 fm to test the 
sensitivity to the lattice artifacts. 

To repeat, a comparison of a dimensionful quantity between quenched and full QCD can 
only be made subject to a choice of the quantity used to determine the length scale in the 
quenched theory. For example, if we choose the string tension, then we might compare the 
wave functions at the origin in units of the string tension, or \l/(0)<7~ 3//4 . 

We solved Schroedinger's equation for valence quark masses of 150, 300, 450 and 750 
MeV. These can be regarded as constituent quark masses for, respectively, a light-light 
meson reduced mass (this is a stretch), a heavy-light meson, a heavy-strange meson, and 
charmonium reduced mass. Figure |7] shows the ground state wave function at the origin 
in units of the string tension for these four masses as a function of the sea quark mass, 
parameterized by (m^/rrip) 2 . As expected, the larger force at short distance results in a 
larger wave function at the origin in full QCD. This is consistent with indications that the 
inclusion of dynamical quarks increases decay constants (fs and Jb s ) |2^-|2B[ and hyperfine 
splittings in quarkonia P J2^J25[1 . This figure also shows ^(O) from the quenched a « 0.10 fm 
run. As expected, the effect of decreasing the lattice spacing is larger for the heavier quark 
masses. As can be seen from Fig. |5], Fig. [7] would look quite different if we chose ro or r± 
to set the scale. For example, at m q = 150 MeV. the wave function in units of ro increases 
by less than one percent as the sea quark mass is lowered to 0.02/a, while in units of a it 
increases by almost five percent. 

An alternative approach is to use the wave function for light quarks to set the scale. This 
is similar to using to fix the lattice spacing, as was done for example in the calculations of 
fs in Ref. |21||. Then we plot the ratios of wave functions for different valence quark masses 



as functions of the sea quark masses in Fig. Again, the dependence on sea quark mass is 
clearly visible, but the size of the effect is much smaller than when the string tension sets 
the scale. For example, ^(O)^ -3 / 4 at my = 300 MeV increases by 6.3 ± 0.5% with three 
degenerate flavors of sea quark with (m^/m^) 2 = 0.46, or about the strange quark mass. 
However, the ratio of wave functions ^300(0) /^isoCO) increases by only 2.2 ± 0.2%. 
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FIG. 7. Wave function at the origin in units of the string tension for valence quark masses of 
150, 300, 450 and 750 MeV. The horizontal scale is the sea quark mass, parameterized by (m^/mp) 2 , 
placing the quenched calculation at the right side. For the two smallest sea quark masses (squares), 
there are two light flavors of sea quark and one heavier flavor at about the strange quark mass. 
The cross indicates the result for the quenched sample with a « 0.10 fm. The vertical scale for 
each graph covers a range of twenty percent. 
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FIG. 8. Ratios of wave function at the origin for different valence quark masses. The horizontal 
axis and symbols are the same as in Fig. The three panels show x I / 3oo(0)/ v I / i5o(0), v I / 45o(0)/ , I'i5o(0) 
and v I / 45o(0)/ l I / 3oo(0), respectively. These can be very loosely interpreted as fs/f-K, fB s /fn and 
Ib s I Ib- Note that in this figure the three panels cover different fractional ranges. 
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